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Structure of Turbulent Shear Flows 

By A.K.M.F. HUSSAIN 1 , J. JEONG 1 AND J. KIM 2 


Introduction 

The accuracy and spatial resolution of numerically simulated databases for tur- 
bulent shear flows, like those generated at NAS A/ Ames, far exceed those typically 
available in laboratory experiments. While there are limitations of the simula- 
tions, in particular regarding low Reynolds number (a technological constraint) 
and limited duration of flow that can be computed (an economic constraint), the 
simulations have also the advantage that they enable the scientists to “measure” 
quantities (such as enstrophy, pressure, dissipation, helicity, and pressure-strain rate 
correlation) which are virtually impossible to measure accurately in the laboratory 
(see Hussain, 1983, 1986). Simulations provide quantitative measures of flow fields 
free from the effects of probe interference and from errors introduced by invoking 
Taylor’s hypothesis. The simulations thus accord heretofore unavailable unique op- 
portunities for research into the structure of turbulence; such was our goal and the 
centerpiece of our effort during the 1987 summer school of Stanford-NASA CTR. 
Our research covered three different topics: 

(1) Eduction of coherent structures, 

(2) Measurement of propagation velocities of perturbations (such as velocity, pres- 
sure, and vorticity) in turbulent shear flows, and 

(3) Direct evaluation of the Taylor hypothesis. 

In the following, we summarize our activities during the summer school in each 
of these three categories. However, most of our effort was devoted to item (1), 
which will occupy the bulk of this report. Recognizing that very little time was 
available to either complete the post-processing or to even adequately digest the 
results obtained, we venture to point out a few apparent interesting observations and 
surprises, make some tentative conclusions and suggest specific areas of continuing 
collaboration between NASA-Ames and University of Houston. 

A. Eduction of coherent structures 

Coherent structures, an embodiment of our search for order in disorder, has 
been the focus of much of the research in turbulent shear flows in the past two 
decades. The overwhelming majority of coherent structure studies has been based 
on flow visualization, which is not only qualitative but can even be grossly mislead- 
ing (Hussain, 1986). We need hard quantitative data regarding the distributions of 
properties over the spatial extent of the structures and the dynamical roles of these 
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structures. These structures being defined in terms of coherent vorticity (Hussain, 
1980), the requirement of vorticity measurement presents a severe constraint. A 
first simplification of this problem was to measure one component of vorticity only 
and in one plane only — that in the azimuthal plane. The initial approach was to 
perform phase-locked measurements triggered on the periodicity of the flow event 
or of a forcing signal (Cantwell &; Coles, 1983; Hussain & Zaman, 1980; Hussain, 
Kleis & Sokolov, 1980). Note that structures in natural (unexcited) shear flows 
undergo jiiiers ox two kinds: initiation and evolutionary jitters. Via excitation, ini- 
tiation of structure formation can be controlled, but the evolutionary jitter remains 
uncontrolled. Thus eduction, even in periodically induced structures, must use a 
local trigger at the measurement location instead of an upstream trigger— a major 
drawback of Cantwell & Coles’ data. 

The next step was to develop a scheme for eduction of coherent structures in an 
unexcited flow by using a local footprint of passing structures (Zaman & Hussain, 
1984). In parallel, efforts were devoted to develop an algorithm that can educe 
structures in any fully developed turbulent flow without requiring any trigger sig- 
nal. This resulted in a scheme which utilizes the measurement signal itself. Such 
a technique (to be explained below) was first developed for eduction of structures 
in the fully developed region of an axisymmetric jet (Tso, 1983) and was further 
refined to educe structures in a turbulent plane wake (Hayakawa & Hussain, 1985). 
This generic and robust scheme, which can be applied to data in any transitional 
and turbulent shear flow, was also successfully applied to numerical simulation and 
experimental data in a plane mixing layer (Metcalfe et al., 1985). The close agree- 
ment between structures educed from numerical and experimental data is strong 
evidence of the robustness of the scheme as well as a validation of the simulation. 

Our goal in this area continues to be: (i) eduction of coherent structure topology 
in different turbulent flows, (ii) understanding of turbulence phenomena in terms of 
entrainment, mixing, production, and dissipation, and (iii) the dynamical role and 
significance of coherent structures in various turbulent shear flows. This provides the 
motivation for eduction of coherent structures from direct simulations of turbulent 
flows. While the coherent structure topology and dynamics in the wall region of 
a turbulent boundary layer would not be expected to be different from those in 
the wall region of a channel, the outer layer structures must be noticeably different 
between the two flows. Hence the motivation for studying coherent structures in a 
flat plate boundary layer and a channel flow. These new data would establish how 
coherent structures in the wall-bounded flows are different from coherent structures 
in free shear flows, which have been studied extensively at University of Houston. 
Having decided to especially focus on the coherent structures near the wall, which 
are responsible for the most interesting events in wall-bounded shear flows, it was 
necessary to separate the effects of the wall from those of shear. Hence the goal of 
eduction of coherent structures in homogeneous shear flows. 

The simulation data in a low Reynolds number, fully turbulent channel flow (Kim, 
Moin & Moser, 1987), in a flat plate boundary layer (Spalart, 1987), and in two 
homogeneous shear flows (Rogers & Moin, 1987; Lee & Reynolds, 1985) are used to 
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educe coherent structures in these flows. 

Eduction methodology: 

Even though the eduction procedure, in principle, is based on the 3D vorticity 
information in a flow field, experimental limitation forced us to implement the 
procedure on one component of vorticity, namely the spanwise vorticity, u> z . In order 
to permit meaningful comparison with experimental data, our earlier eduction from 
numerical simulation was also based on only the w z field of the simulation (Metcalfe 
et al., 1985). 

The eduction procedure is briefly outlined here. Readers can find further details in 
our earlier papers. The eduction steps are as follows: (i) select a small y - range where 
coherent structures are to be studied and obtain u > 2 (x,y) data from simulations in 
various 2 -planes, but centered at the middle of this y- range (x, y, z are longitudinal, 
transverse, and spanwise coordinates, respectively); (ii) smooth the a > 2 (x,y) data 
using a zero-phase shift filter; (iii) detect structures which are strong (i.e., peak 
vorticity above a threshold w tl ) and of a sufficiently large size; (iv) look for z- 
symmetry of these vorticity concentrations; (v) accept u> 2 data from a plane of 
symmetry (being careful not to accept two planes of the same structure). 

The phase average of all realizations containing similar structures yields a coher- 
ent structure. Note that a coherent structure is a stochastic quantity and may not 
be observed instantaneously. The center (marked by the peak value) of smoothed 
w z contours is only a first guess, and is a useful reference for initial alignment of 
various realizations, but may have nothing to do with the true center, which must 
be determined from the ensemble averaging after proper alignment. 

The eduction continues as follows: (vi) align w z with respect to the peaks of 
smoothed contours and obtain the ensemble average — this is the zeroth iteration 
ensemble average D 0 ; (vii) obtain cross correlation 0 of each realization with 
the ensemble average; (viii) shift the center of each realization by the location in 
(x,y) for peak correlation; (ix) reject realizations requiring excessive shifts; (x) 
obtain ensemble average of the realizations by aligning them with respect to the 
new centers — this is the first iteration ensemble average ; (xi) reject realizations that 
produce excessively low correlation peaks; (xii) continue the iteration until all shifts 
required for alignment fall below a size and all correlation peaks are above a set 
level; (xiii) identify the locations of revised centers of finally accepted realizations 
in the unsmoothed data records and obtain ensemble average after aligning with 
respect to these centers — this final ensemble average is the coherent structure; (xiv) 
the departure of each unsmoothed realization from the ensemble average denotes 
incoherent turbulence. 

In the eduction, we set u) tl in terms of the local maximum of mean u? 2 (y). In case 
of the channel flow, the structures educed were those centered at y + — 125, 50, 30, 
and 15. 

The educed velocity vector patterns in the (x, y)-plane for the four locations are 
shown in Figs, l(a-d). All the four contours show some similarities, characterized by 
saddles and centers. We had anticipated the structure at y + = 125 (note y + = 180 at 
the centerline of the channel) to resemble to some extent the mixing layer structure 




FIGURE 1. Educed velocity vector patterns in (x,y)-plane: (a) Centered at 
y+ = 15; (b) y + = 30; (c) y+ = 50; (d) y+ = 125. Note that the y-scales are 





Structure of Turbulent Shear Flows 


277 


studied extensively before (see Hussain & Zaman, 1985; Metcalfe et al., 1985), but 
the near wall structures appear to be different presumably due to the effects of the 
wall and high shear. For this reason we present detailed data for structures centered 
at both y + = 125 and y + = 15. The results for y + = 15 are presented in Figs. 2(a- 
m) as contours of coherent vorticity < u> 2 >; coherent longitudinal and transverse 
velocities < u > and < v >; coherent pressure < p >; incoherent turbulent kinetic 
energy < q 2 r >; incoherent Reynolds stress -< u r v r >; incoherent pressure < p 2 r >; 
coherent strain rate < S > = <J)u/dy > + < dv/dx >; coherent shear production 
< Ps >= ~ < U r v r >< S >; coherent noripal production < P N >= - < u 2 r > d < 
u > /dx— < v 2 > d < v > /dy\ total production < Pt > = < Ps > + < Pn >; 
coherent pressure works < p >< S u >, < p >< S 22 >. Note that the y- scale 
is expanded relative to the x-scale in order to reveal the details of the flow field 
and that two vorticity contours are duplicated in all figures to provide a common 
reference. The corresponding contours for y+ = 125 in the channel flow are shown 
in Figs. 3(a-m). 

We have started educing structures at four transverse locations in the Spalart’s 
boundary layer as well as one location each in the homogeneous shear flows. For 
the purposes of this preliminary report, we limit our discussions primarily to the 
case of channel structures centered at y + = 15. 

In the frame of the advected coherent structure, the flow above it moves down- 
stream and the flow below it moves upstream. The structure advection velocity 
being about 60% of the centerline velocity, the downstream stagnation point is 
closer to it than the upstream one. One would thus expect the normal production 
at the front to be higher than at. the back, as is indeed the case (Fig. 3j). Note 
that nearer to the center of the channel, the shear is weak and more nearly uniform 
across the structure. That is why the two saddles are nearly equidistant from the 
center, and the normal productions are equal on both the front and the back. 

The coherent pressure contour extends in the transverse direction considerably 
beyond the structure boundary as denoted by coherent vorticity. This is to be 
expected from the fact that pressure is an integral property, being the solution of a 
Poisson equation with the source term due to gradients of the velocity field. Note 
that the minimum of coherent pressure is at around the structure center, but does 
not exactly coincide with it. 

One striking feature is the fact that longitudinal pressure work is mostly negative 
(hence of the right sign) at y + = 125 but mostly positive at y+ = 15 (hence of 
the wrong sign). That is, the pressure work transfers kinetic energy away from 
the longitudinal component in the outer layer as is commonly presumed on the 
basis of time average kinetic energy balance (see Tennekes & Lumley, 1972), and it 
transfers kinetic energy into the longitudinal component (contrary to expectation) 
near the wall. This is consistent with the result of Moin & Kim (1982), in which 
they attributed it to the “splatting” motions of large eddies in the near- wall region. 

Since the shear layer below the advecting structure (near the wall) has a higher 
velocity gradient than the one above the structure (toward the centerline), one 
would expect the shear production < P s > to be higher on the left- than on the 
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right-hand side of the structure. Data show the two regions to have comparable 
levels of < P s >, suggesting that the incoherent Reynolds stress < u T v r > is higher 
at the top right-hand side of the structure. Note that both < Pjy > and < Ps > are 
coordinate dependent and it is their sum which is invariant under rotation (Hussain, 
1983). The total production < Pt > is higher at the top on the right- hand side. 

Concluding remarks 

Tb*> success of the eduction scheme in cAiiaciing coherent structure details in 
the fully turbulent channel flow is demonstrated. In the plane of ^-symmetry the 
coherent structure characteristics are quite similar to those in the free mixing layer — 
more so in the outer regions than near the wall. The topology consists of saddles 
and centers, the saddle region being the location of maximum incoherent Reynolds 
stress -< u r v r >, and maximum shear production. One interesting difference from 
the mixing layer case is that the center in the wall-bounded case is not necessarily 
characterized by a high level of incoherent turbulence intensity. 

While we need to devote more time to compare our channel data with those in 
jets, wakes and shear layers, some of the measures have not been obtained yet in 
the free shear flows: for example, contours of dissipation, pressure work, etc. These 
measurements and simulations are planned in the future. One noticeable difference 
between the educed structures at = 125 and y + = 15 is that the structure in the 
former case is more rounded while in the latter case it is much more sheared. In 
order to determine the role of shear while eliminating any effect of the wall, it would 
be worthwhile to educe these structures in homogeneous shear flows with different 
shear rates. This work is in progress now. 

While we do not expect any noticeable difference in the near- wall coherent struc- 
ture characteristics between channel flow and flat plate boundary layer, the outer 
layer structures can be quite different in the two flows on two accounts. In the 
boundary layer, the outer structures should all be of the same sign and be bounded 
by irrotational (nonturbulent) fluid which they entrain. In the case of the channel, 
the outer region should consist of structures originating from both walls and can 
be significantly different as a result of interaction of these structures of oppositely- 
signed circulations. 

When all these data are completed, we will be able to identify the similarities 
and differences among structure topologies in the channel flow, the boundary layer 
and homogeneous shear flow. We then hope to be able to comment specifically on 
the role of the nonturbulent freestream, on the role of shear, and on the role of the 
wall. 


Future extensions 

We propose to extend this work to include the following: 

(1) Starting from the plane of symmetry, march ahead on either side to track the 
coherent vortical structure and educe structure details in local planes and then 
reconstruct the three-dimensional structure. 

(2) Educe structures with circulation opposite to the direction of mean shear. This 
is important in the outer region of the channel where structures migrate from 
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Figure 2. Educed turbulence structures centered at y+ = 15. Note that the 
y-scales are stretched to reveal the details of the structures, (a) < u> 2 >; (b) 

< u >; (c) < r >; (d) < p >; (e) < q 2 r >; (f) -< Ur v r > ; ( g ) < > ; (h) 

< S > = < du/partialy > + < dv/dx >; (i) < P s >= - < Ur v r >< S >; (j) 

< P N >= - < ul > d < u > jdx- <vj> d<v> jdy\ (k) < P T >=< P s > + < 

Pn >; (1) < p >< Su >; (m) < p >< S 2 2 >• 


ORIGINAL PAGE IS 
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FIGURE 2. ... continued. 


the opposite wall. 

(3) Complete eduction on the basis of total vorticity |u;| rather than the spanwise 
vorticity u> 2 . 

(4) Address clearly the mechanisms of entrainment, mixing, production, intercom- 
ponent transport, enstrophy cascade and dissipation associated with the educed 
three-dimensional coherent structure. 

(5) Educe other significant coherent structures and evaluate the dynamical role and 
significance of coherent structures in the three flows studied. 

(6) Study the evolution of a turbulent hairpin vortex in a laminar boundary layer 
and a laminar channel flow with a velocity profile matching that of a turbulent 
flow. Study the same in a flow with laminar and turbulent homogeneous shear. 

(7) Study the dynamics and evolution of an an artificially induced bursting coherent 
structure in fully turbulent channel flow and flat plate boundary layer. 

B. Propagation velocities 

The direct simulation data provide an excellent opportunity for determination of 
propagation velocities of pressure, velocity and vorticity perturbations in turbulent 
shear flows. In the case of the turbulent boundary layer or channel flow, in partic- 
ular, varying values of propagation velocities have been reported in the literature. 
In addition to being of fundamental interest, the propagation velocity is of direct 
concern in understanding coherent structure topology and dynamics. 

If instantaneous fields of simulation data are considered at two instants <j and <2 
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Figure 3. Educed turbulence structures centered at y+ = 125. Note that the 
y- scales are stretched to reveal the details of the structures, (a) < u> 2 >; (b) 

< « >; ( c ) < v >; (d) < p >; (e) < qj >; (f) -< u r v r >; (g) < p 2 r >; (h) 

< S > = < du/partialy > + < dv/dx >; (i) < P s >= - < Ur v r X S >; (j) 

< P N >= - < u 2 r > d < u > jdx- < v 2 > d < v > /dy; (k) < P T > = < P s > + < 
P N >; (1) <p>< S u >; (m) < p >< 5 22 >. 
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Figure 3. ... continued. 


and they are cross-correlated, then from the locations ( 6x m ,Sy m ) of cross- correla- 
tion peak, one can obtain streamwise and transverse propagation velocities U c , V c 


as follows: 


Uc 


tern 

(* 2-U) 


v c = 


hjn 

(t 2-U) 


As a first step, we have determined values of U c from streamwise correlations of 
velocity («i), pressure (p), and vorticity (u>i) fields. Fig. 4a shows the profiles of 
U c for Ui and p as functions of y, and Fig. 5a shows the profiles of U c for vorticity 
(u>i). Figs. 4b and 5b show the same profiles as functions of y + . 

It is surprising how closely the convection velocities for velocity and vorticity 
perturbations agree with each other. Moreover, these profiles closely agree with 
the mean velocity profile, being only slightly lower than mean velocity in the outer 
region, but being higher than the mean velocity near the wall. There are clear 
differences between the data presented here and those reported in literature. The 
convection velocity of pressure in the wall region is consistently higher than those of 
velocity and vorticity, indicating somewhat the elliptic nature of the pressure field. 


Future extensions 

(1) To determine V c in addition to U c . When correlation in (x and y) is used, 
the current value of U c may be somewhat different from that found from the 
x-correlation alone. 
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(2) Determine U c and V c for positive and negative w z separately. 

(3) Determine U c and V c values for all variables in the two homogeneous shear flows. 

C. Taylor hypothesis 

The Taylor hypothesis of frozen turbulence is a common assumption in experi- 
mental turbulence for inferring spatial structure of turbulence from temporal data 
obtained by a stationary probe. It is also commonly invoked for inferring wave 
number spectrum from measured frequency spectrum as well as in measurements of 
dissipation and higher-order moments. While there have been extensive discussions 
of the limitation of the hypothesis, which were obviously known to Taylor himself, 
no direct test of this hypothesis has been possible yet. The direct numerical sim- 
ulation databases allow us to make a thorough evaluation of this hypothesis for 
different flow variables as a function of shear rate. 

The hypothesis assumes that 

D T* = l§i + Vt£)*=0 

where Ut is the Taylor advection velocity and <j> can be any of the variables 
{w,,p,u> t , c}. The value of U T used in the literature has varied (Zaman & Hus- 
sain, 1981): 

Ut = U (y), mean velocity 

= u(x,y,z), local velocity 

= U c (y), propagation velocity 

= Uf(x, y, z), filtered velocity 

We have used the first, second and fourth in the 1987 CTR session. It is seen that 
both mean and rms values of Dr<i> are small in the outer layer, except very close to 
the wall. Comparing instantaneous contours of D T <f> with the smoothed contours of 
u i , p, it is found that the departure from the hypothesis is not directly associated 
with large-scale structures. 

Future extensions 

(1) Compute vorticity fields using Taylor hypothesis and velocity fields, and com- 
pare with true vorticity fields. 

(2) Determine which choice of Ut produces the minimum error in the use of the 
Taylor hypothesis. 

(3) For situations when there are large values of D T <t>, evaluate the various terms 
in the convective balance equations of <f>. Determine which term (and hence 
which mechanism) produces the maximum contribution to the departure from 

D t 4> — 0. 
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